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ABSTRACT 

We derive Galactic continuum spectra from 5-96 cm~^ from COBE/FIRAS 
observations. The spectra are dominated by warm dust emission, which may 
be fit with a single temperature in the range 16-21 K (for emissivity) along 
each line of sight. Dust heated by the attenuated radiation field in molecular 
clouds gives rise to intermediate temperature (10-14 K) emission in the inner 
Galaxy only. A widespread, very cold component (4-7 K) with optical depth 
that is spatially correlated with the warm component is also detected. The 
cold component is unlikely to be due to very cold dust shielded from starlight, 
because it is present at high latitude. We consider hypotheses that the cold 
component is due to enhanced submillimeter emissivity of the dust that gives 
rise to the warm component, or that it may be due to very small, large, or fractal 
particles. Lack of substantial power above the emission from warm dust places 
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strong constraints on the amount of cold gas in the Galaxy. The microwave sky 
brightness cItic to interstellar dust is dominated by the cold component, and its 
angular variation could limit our ability to discern primordial fluctuations in the 
cosmic microwave background radiation. 

Subject headings: Infrared: Interstellar: Continuum — ISM: Dust, Extinction 
— Radio Continuum: Interstellar 



1. Introduction 

Absolute spectrophotometry on large angular scales in the far-infrared (A ^ lOO/xm) 
and submillimeter (A ^ 1mm) bands is not possible using ground-based instruments, 
because of the large, time-varying attenuation by the Earth's atmosphere. Even using 
balloon-based instruments, observations require frequent chopping to 'blank' regions, so 
that the accurate photometry is available only on scales comparable to the reference 
beam separation. Portions of the galactic plane were surveyed by balloon instruments in 
wavebands between 150 and SOO/um, revealing the bright ridge of dust emission from H II 
regions and molecular and atomic clouds centered on the galactic plane (cf. Hauscr ct al. 
1984 and references therein). The Infrared Astronomical Satellite (IRAS) made the first 
all-sky survey at 100/i.m wavelength (Neugebauer et al. 1984). The lOO/xm emission at high 
galactic latitude is strongly correlated to the distribution of atomic gas as traced by the 
21-cm fine of H I (Boulanger & Perault 1987), and can be explained by emission from dust 
grains, mixed with the interstellar gas and heated by the interstellar radiation field (Draine 
& Lee 1984). Interstellar grains of various sizes contribute different parts of the infrared 
emission. In one recent model, the emission is broken into three components, due to 'big' 
grains, 'very small' grains, and polycyclic aromatic hydrocarbons (Desert, Boulanger, & 
Puget 1990). The far-infrared emission is dominated by 'big' grains, which also dominate 
the mass and total luminosity of interstellar dust. 

In this paper we describe all-sky observations by the Far Infrared Absolute 
Spectrophotometer (FIRAS) aboard the Cosmic Background Explorer (COBE), which 
provide the first complete coverage of the wavelength range from 4.5mm through 104/im. 
Preliminary results on the Galactic continuum and spectral lines in this wavelength range 
were presented by Wright et al. (1991), and the angular variation of the brightnesses of 
the spectral fines was presented by Bennett et al. (1994). Here we focus on the galactic 
continuum spectrum and its variation both in the galactic plane and at high galactic 
latitude. 
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2. Observations 

2.1. FIRAS Description 

The FIRAS is one of the three Cosmic Background Explorer (COBE) instruments. 
During its observing Ufetime (November 1989 to September 1990), 95% of the sky was 
covered by scans at a solar elongation of 91-94°. The effective beam has a FWHM of 7° 
and drops sharply to -30 db by 10° off-axis, as confirmed by observations of the Moon. 
The spectrometer is a Michelson interferometer with bolometric detectors sensitive to 1-96 
cm~^ radiation. (Throughout this paper, we use the cm~^ scale for frequencies, u. To 
convert to GHz, multiply by 30; to convert to wavelength in //m, use 10'^ /u.) The incoming 
photons are separated into high- and low-frequency channels by a dichroic filter. The 
spectrometer was operated in several modes by varying the stroke length and scan rate. 
For this paper the data from the low and high frequency detectors were combined, so that 
the spectra cover frequencies 2.2-96 cm^^ (or wavelengths 4.5mm-104/im), with a spectral 
resolution of 0.57 cm~^ over the entire range. (See Fixsen et al. 1994 for technical details). 
The observations were absolutely cahbrated by placing a beam-fiUing blackbody in the 
aperture during calibration periods (once per month during most of the mission). Relative 
calibration was maintained by always measuring the difference between the sky signal and 
that of an internal calibrator. 



2.2. Data Prepeiration 

The monopole and dipole components of the cosmic background radiation were removed 
from each spectrum by subtracting a blackbody with a temperature Tqer = + ^dCOsG, 
where Tq = 2.726 K, = 3.343 mK (Mather et al. 1993), and is the angle between 
the line of sight and the orientation of the cosmic dipole (/ ~ 264.4°, h — 48.4°) (Kogut 
et al. 1993). The resulting spectra are dominated by the continuum emission of galactic 
dust, spectral lines of galactic gas, the (as-yet undetected) cosmic infrared background, 
and contributions from the zodiacal light at high frequencies and galactic synchrotron and 
free-free emission at very low frequencies. 

The zodiacal fight was modeled using observations by the COBE Diffuse Infrared 
Background Experiment (DIRBE). The DIRBE observations at 5 wavelengths between 4.9 
and 100//m were used to fit a parameterized model of the distribution and temperature 
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of interplanetary dust. The solar elongation of the DIRBE line of sight varied from 
64° to 124° during each spin about the COBE rotation axis, introducing a strong time 
dependence to the sky brightness. The parameterized model was optimized to match the 
temporal variation of the sky brightness; therefore, no assumptions were made about the 
time-invariant sky brightness outside the Solar System. The parameterized zodiacal light 
model was evaluated for each FIRAS line of sight, and for a frequency of 70 cm~^, to 
produce a spatial template, Z{l,b). This spatial template was then used to determine 
the zodiacal light spectrum, by assuming the brightness can be decomposed into cosmic. 
Galactic, and zodiacal components: 

b) = B,{T,,r) + 9uG{l, h) + z,Z{l, b), (1) 

where the cosmic microwave background is assumed to be the only isotropic background. 
The decomposition of the FIRAS spectrum for the region (/ = 20 ± 5°, 6 = 20 ± 5°) is shown 
in Figure This technique is similar to that used by Wright et al. (1991) to determine the 
spectrum of the Galaxy, gi^; for this paper, we excluded the galactic plane {\b\ < 30°) and 
added the zodiacal term. The resulting zodiacal light spectrum is adequately represented 
by a power-law, ^(z/) oc u^, over the FIRAS bandpass. We removed the zodiacal light 
from each FIRAS spectrum by subtracting (z//70)^Z(/, 6). If the zodiacal spatial template 
is accurate to better than 10%, then the zodiacal light subtraction introduces negligible 
uncertainty except at the very highest frequencies and at the galactic pole. At high galactic 
latitude {\b\ > 60°), the zodiacal light contributes 60% of the total emission at 100 cm~^, 
25% at 70 cm-\ and 10% at 42 cm'^ 

The contributions of Galactic synchrotron and free-free emission to the FIRAS 
frequencies were constrained using observations by the COBE Differential Microwave 
Radiometer (DMR). The DMR observations at three frequencies, 31, 53, & 90 GHz (1.0, 
1.8, & 3.0 cm^^), were extrapolated to higher frequencies assuming the spectrum follows a 
power-law, S,y oc u'^. For synchrotron emission, the spectrum follows /? ~ 0.75 at lower 
frequency (cf. de Bernardis, Masi, & Vittorio 1991), and then the spectral index steepens 
to /3 ~ 1 above 10 GHz (Banday & Wolfendale 1991). For free-free emission, the spectral 
index is /3 = 0.1. In order to place an upper limit to the free- free and synchrotron emission 
at higher frequencies, we used the shallower /5 = 0.1 to scale the total brightness at 53 
GHz. This is very conservative, because a substantial fraction of the 53 GHz emission is 
due to synchrotron and dust emission. A cosecant of Galactic latitude was fitted to the 
53 GHz map, and an isotropic component with a brightness equal to the cosecant slope 
(12/xK) was added to the map to replace the isotropic component (cf. Bennett et al. 1992). 
The combined synchrotron and free-free emission account for < 10% of the emission at 
frequencies higher than 5 cm~^. In what follows, we excluded frequencies below 5 cm~^ in 
order to avoid any possibility of contamination by free-free emission. 
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There were 9 spectral lines detected by FIRAS in the average spectrum of the Galaxy 
(Wright et al. 1991, Bennett et al. 1994), and several of these are evident in each 
spectrum at low galactic latitude. These lines were avoided for the purposes of this paper 
by eliminating the appropriate frequencies. All of the lines are unresolved, so only two 
frequency bins contain each line emission; however, some ringing in the apodization function 
was evident for the very bright lines. We removed spectral frequencies within 1.5 cm~^ of 
the bright C II (63.395 cm'^) and N II (48.72 and 82.04 cm'^) lines, and within 1 cm'^ of 
the C I (16.419 and 27.00 cm'^) and CO (7.681, 11.531, 15.375, and 19.237 cm'^) lines. 

Coaddition was necessary in order to accurately measure the galactic spectrum at 
high latitudes. The spectra were combined into 120 longitude bins in the galactic plane 
and 26 high-latitude regions. The bins were chosen to be larger at higher latitude, where 
the signal is relatively weak. For the galactic plane spectra, the bins are evenly spaced 
bins at 3° longitude intervals and include data with |6| < 3°. The high latitude bins 
consist of the north and south galactic polar caps, with \b\ > 60°, high-latitude zones 
with 60° > |6| > 30° in 90°-wide longitude bins, and intermediate- latitude zones with 
30° > |6| > 10° in 45°-wide longitude bins. Two sample spectra, toward (/ = 45°, 6 = 0°) 
and (180° > / > 90°, —30° > b > —60°) are shown in Figures ||(a) and 0(a). The error bars 
shown here and used in the spectral analysis are the calibration uncertainties in the mean 
brightness over each spatial region, and they do not reflect the intrinsic variation of the 
brightness within the region. 



2.3. Spectral Analysis 



The spectrum for each region was fitted by one or more modified blackbodies of the 

form 

I, = Toe,B,{T), (2) 

where tq is the optical depth at = 30 cm^^ (Aq = 333/im), and ei, is the emissivity 
normalized to unity at z/q. The simplest fit for each spectrum was a single component 
with a power-law emissivity, e,y = (u/uq)"'. In the low-frequency limit, the Kramers- Kronig 
theorem can be used to show that a is an even positive integer (Wright 1993), although it 
is not yet known below which frequency this limit applies. The simplest fit we considered 
was for a u"^ emissivity, so that there are only 2 free parameters (T and r). As examples, 
the residual spectra following two such fits are shown in Figures 0(b) and 0(b). The 
residuals clearly indicate the presence of excess emission at low frequencies in nearly all of 
the spectra, indicating a very cold emission component that is widespread. We repeated the 
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spectral fits allowing a to vary, which substantially improved the of most fits. Examples 
of the residuals that resulted in this case are shown in Figures 0(c) and 0(c). 

In the inner Galaxy, where the signal to noise is high, poor fits were obtained with a 
power-law emissivity. The galactic center spectrum at frequencies below 50 cm~^ is best 
fit by a = 1.98 ± 0.01, and the other spectra require lower a. There are good reasons to 
expect that the value of a derived from fits using single temperatures is lower than the true 
value: the distribution of temperatures along the line of sight should broaden the spectrum, 
and lower our derived a. The highest observed value of a is most likely to approach the 
true value, and we will therefore adopt a = 2 at low frequencies. At higher frequencies, the 
emissivity increases less steeply than . A broken power-law model that can match the 
emissivity rolloff at high frequencies has a = 1.5 for frequencies above 70 cm^^; however, 
the broken power-law introduces an unobserved cusp so that is large for these fits. A 
smooth model for the emissivity variation that can match the high-frequency rolloff and 
improve the of fits is 

= r, (3) 

with z/i ~ 50 cm~^ and /3 ~ 6. These values of ui and the exponents in equation |] were 
chosen by experiment, and are not the result of a formal fit; they improve for the / = 15° 
fit by a factor of 2. The high-frequency rolloff is needed only in the inner galaxy, but is not 
excluded in the high-latitude spectra, where the signal-to-noise is low. At high latitude, a 
simple power-law (z/^) is adequate. It cannot be determined with the present data whether 
the high-frequency rolloff in the emissivity is due to different grain properties in the inner 
galaxy, or the distribution of grain temperatures along the line of sight. In what follows, we 
use equation || for the galactic plane and a simple z/^ for the high-latitude spectra. 



3. Results 

3.1. Significance of the spectral models 

The FIRAS spectra of galactic emission can be characterized by three components, 
which we call 'warm' (16 < Ti < 21 K), 'cold' (4 < T2 < 7 K), and 'intermediate' 
(10 < T3 < 13 K). There is no overlap between the temperatures of the three components, 
and the ranges are those of the fit results. Each spectrum was fitted by a single-component, 
a = 2 model (T, r varying); a single- component, free a (T, r, a varying); and a 
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two- component, a = 2 model (Ti, ri, T2, T2 varying), with the goal of finding the simplest 
model that could fit the most spectra. The significance of a given model was judged by 
its improvement in the reduced x^) according to the F-test. The probability that the 
improvement in due to including an extra parameter was not due to chance, P{> F), 
was calculated. In every case, the single-component, free-a model and the two-component 
models were significantly better than the single-component, a = 2 model, and in most cases, 
the two-component model was significantly better than either single-component model. In 
Table 0, the results of the two-component model for the galactic plane spectra are listed, 
together with the probability, P{> Fac), that the two-component model is significantly 
better than the single- component, a = 2 model, and the probability, P{> Fhc), that the 
two-component model is significantly better than the single-component, free a model. For 
some spectra, it was not possible to distinguish between the single-component, free a and 
two-component models at the 90% confidence level. In the inner Galaxy, both the one- 
and two-component fits are poor because of the intermediate-temperature emission. A 
three-component fit improves iii the inner galaxy, but cannot be fitted to the outer 
galaxy or high-latitude spectra because the third component cannot be separated from the 
other two. We therefore consider the two-component model for the analysis of the galactic 
plane spectra. 

The results of some of the high-latitude fits are shown in Tables b, and c for 
the various spectral models. The excess submillimeter emission above a u"^ modified 
blackbody is present with > 95% confidence in all spectra. The distinction between the 
single-component, free a model and the two-component model was not possible for many 
high-latitude spectra. The results of the free-a fit indicate statistically significant variations 
in a, so that the emissivity law varies substantially from place to place. Based on the fact 
that the two-component model was more successful in the galactic plane, it is likely that 
the modest amount of cold emission, even at |6| > 60°, is also due to a second component. 



3.2. Effect of a cosmic infrared background 

The cosmic infrared background radiation (CIBR) was assumed to be small in this 
analysis, although the observational upper limits to its brightness still allow it to contribute 
a significant fraction of the sky brightness. Models of the integrated light of galaxies 
(Franceschini et al. 1994) predict CIBR brightnesses that are 5-25% of the 10-100 cm~^ 
sky brightness (after zodiacal light subtraction) at the galactic poles. In order to determine 
whether the presence of a CIBR would have a significant effect on our conclusions, we 
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subtracted the brightest CIBR from Franceschini et ah (1994) from the FIRAS spectra and 
repeated the spectral fits. The optical depths decreased, as expected because the sky is 
fainter, but the temperatures of the warm and very cold components are not significantly 
affected. We therefore feel confident in neglecting the CIBR, noting only that the dust 
optical depth at the highest latitudes is overestimated. 



Two effects will produce a range of temperatures in each spectrum. First, each particle 
size has its own temperature, determined by the balance between its far-infrared and visible 
absorption efficiencies. Second, different regions of the Galaxy have different radiation 
fields, which heat the dust to different temperatures. The issue of disentangling these 
effects is beyond our present scope. The radial variation of the interstellar radiation field 
can best be determined by comparison to spectral-line surveys of CO and H I, for which the 
distances can be determined from Galactic rotation. In order to gain some insight into the 
effect of a temperature distribution, we examined the following spectral model: 



which assumes equal probabilities for grains at each temperature in the range To ± AT. 
There is a strong covariance between T and AT, such that opening the width of the 
temperature distribution leads to lower mean temperatures. We have fitted eq. § to 
representative FIRAS spectra for the inner Galaxy, outer Galaxy, and the high-latitude sky. 
For the inner Galaxy, the optimum width to the temperature distribution was AT = 0. 
The fact that eq. |^ does not improve the fits, despite adding an extra degree of freedom, 
leads us to conclude that either the temperature distribution is not important, or it has 
a substantially different functional form than we considered. For the outer Galaxy, eq. ^ 
produced fits of comparable quality (x^) for a wide range of AT, as long as the maximum 
temperature, Tq + AT ^ 19 K. Thus it is possible that comparable column densities of dust 
exist at a wide range of temperatures in the outer Galaxy. The power emitted depends 
approximately on the 5-6*^ power of temperature, so negligible amounts of power are 
emitted by dust much cooler than that at the highest temperatures. At high latitudes, 
eq. eqdis was no better than a single temperature. It is worth noting that changing AT has 
only a small effect on the shape of the spectrum at low frequencies. Thus using a single 
temperature, rather than a continuous distribution, for the warm dust emission has no 
impact on the presence of the 'very cold' component in the spectra. 



3.3. Effect of a temperature distribution 




(4) 
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The small width of the temperature distribution for the inner Galaxy has direct 
implications for the variation of temperature with galactocentric radius and particle size. 
An upper limit to the width of the temperature distribution for the line of sight in the 
Galactic plane toward / = 45° is AT < 2 K at the 95% confidence level. The expected 
AT due to Galactic radial variations may be estimated by assuming the volume emissivity 
scales exponentially with scale length R^, and and the radiation field has scale length i?^,; 
then AT/T ~ RJ6R^. If R, = 3.1 kpc, R^ = 3.5 kpc, and T ~ 16 K, then AT ~ 2 K, 
comparable to our upper limit to AT. The expected AT due to the range of particle sizes 
may be estimated as the rms dispersion in T weighted by the emissivity per particle size; 
we find AT ~ 1 K from this effect. Thus the bulk of the far-infrared emission is fit well 
by a narrow temperature range, and the expected temperature range due to the Galactic 
radial gradient and the size distribution is comparable to the upper limit derived from the 
data. This justifies the assumption of a single temperature along the line of sight as a first 
approximation. 



3.4. Distribution of temperature and optical depth 

The longitude profiles of the optical depths and temperatures of the warm and cold 
emission are shown in Figure ^ It is evident that the temperature of the warm component 
varies strongly with galactic longitude. This variation is expected because lines of sight 
passing closer to the galactic center sample dust subjected to a higher radiation field 
(Mathis, Mezger, & Panagia 1983). The temperature of the cold component does not show 
the same longitude variation; if anything, the temperature of the cold component may be 
lower in the inner Galaxy. 

The high-latitude variations of the temperature and optical depth of the 'warm' and 
'cold' components are shown schematically in Figure |^. In the 10° < |6| < 30° zone, the 
temperature of the warm component generally follows that in the galactic plane, with clear 
exceptions being the relatively hot Orion starforming region at (/ = 210°, 6 = —20°) and 
the relatively cold Taurus dark cloud region at (173°, —14°). Note that the optical depth is 
high in both the Orion and Taurus regions, demonstrating that the temperature and optical 
depth vary independently and can be separately determined. The Orion region contains a 
large OB association, and has a radiation field enhanced by a factor of ^ 2 with respect to 
the general interstellar radiation field, when averaged over scales comparable to the FIRAS 
beam (Wall et al. 1994). The Taurus region has relatively few OB stars, and its infrared 
emission is dominated by relatively dense, cold molecular clouds. At higher latitudes, the 
temperature approaches the polar value of around 18 K, with both poles comparable. 



- 10 - 



The optical depths of the warm and cold components are well correlated. In Figure ^ 
they are plotted against each other. The rank correlation coefficient between ti and T2 is 
0.96 for all spectra, and it is 0.95 for the high-latitude spectra only. The correlation is 
not quite linear, with the average T2/T1 = 7.1 ± 0.4 for high latitudes (low optical depths) 
and 5.2 ± 0.2 for low latitudes (high optical depths). It is remarkable that the cold and 
warm components are so well correlated, and that nearly the same ratio of optical depths 
is obtained for all lines of sight. However, it is also significant that the correlation is not 
perfect, so there is not a single, universal type of dust. 



3.5. Comparison to Previous Observations 

The distinction between the warm and very cold components was made possible by 
the continuous spectral coverage of FIRAS, extending throughout the far-infrared and 
submillimeter bands. Observations in the 100-350/im range are consistent with dust at 
a single temperature — the warm dust. The warm dust temperatures we measure in the 
galactic plane are similar to those obtained from balloon-based surveys (cf. Hauser et al. 
1984), and the high-latitude dust temperature we measure is similar to that (16.2ti j K) 
obtained from broadband (134, 154, & 186/im) observations by a Berkeley/Nagoya rocket 
experiment (Kawada et al. 1994). The very cold emission has been detected before, and 
ascribed either to a submillimeter excess or to an emissivity index a < 2. Combining 
broad-band observations at 5.6, 8.7, 15.8, and 22.5 cm~^ (1800, 1100, 630, and 440 /im) by 
a MIT balloon-borne instrument (Page, Cheng, & Meyer 1990) with observations at 33, 40, 
and 67 cm~^ (300, 250, and 150 /xm) by a Goddard balloon-borne instrument (Hauser et al. 
1984) and observations at 100 and 167 cm~^ (100 and 60 /xm) by IRAS, the spectrum of 
the galactic plane at / = 42° was found to have a submillimeter excess (Page et al. 1990). 
The excess could be characterized by a temperature of 4 K and optical depth 4 x 10~^ at 30 
cm~^. We have compared their observations to the FIRAS spectrum at the same location. 
The excess submillimeter emission from the MIT/Goddard spectrum is larger than that 
seen by FIRAS, but the results are consistent to within the uncertainties due to different 
beam sizes and calibration. Similarly, combining broad-band observations at 6, 9, and 12 
cm~^ (1700, 1100, 830 fim) by a Santa Barbara/Berkeley balloon experiment (Meinhold 
et al. 1993) with IRAS 100/im observations, the best-fitting single-component model for 
the dust emission required a = 1.4 for lines of sight near the star /i Peg (6 = —30.7°). A 
low value of a indicates excess emission with respect to the z/^ emissivity expected at low 
frequencies, and therefore the Santa Barbara/Berkeley results provide more observational 
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evidence for the very cold component in the spectrum of interstellar dust at high galactic 
latitude. 



4. Discussion 

4.1. Warm Dust (16-23 K) 

The temperature of the warm component of the FIRAS spectra agrees with predictions 
for dust with a size distribution consistent with interstellar extinction (Mathis, Rumpl, & 
Nordsieck 1977) and heated by the interstellar radiation field (Mathis, Mezger, & Panagia 
1983; Draine & Anderson 1985). They are also similar to those derived from broadband 
60-100 fim observations by IRAS (Sodroski et al. 1987; Bloeman, Deul, & Thaddeus 1990). 
However, the temperatures inferred from the 60-100 fim ratios do not show the longitude 
variation expected for dust heated by a radiation field that increases toward the galactic 
center. This mystery, which has posed problems for models of the physical properties of 
the dust grains responsible for the emission (Desert, Boulanger, & Puget 1990), is resolved 
now that all-sky surveys at longer wavelengths are available. The DIRBE broadband 
140/240 /im temperature (Sodroski et al. 1994) shows the same longitude variation evident 
in Figure ^. The FIRAS observations demonstrate that the spectrum of interstellar dust 
emission in the 100-300 /im range is due to emission from grains with nearly a single average 
temperature. The excess mid-infrared emission, which dominates shortward of 40 fim 
and contributes strongly at 60 /im, is evidently due to particles undergoing temperature 
fluctuations (Draine & Anderson 1985; Desert et al. 1990). 

The fact that the interstellar dust spectrum between 100-300 /im can be fitted by 
a single modified blackbody suggests that it is produced by large grains in equilibrium 
with the interstellar radiation field. Thus observations longward of 100 /im can be reliably 
used to infer the optical depth of interstellar dust. The narrow temperature range also 
suggests that the emission is not produced by multiple grain types with widely different 
temperatures. For example, a model with bare graphite and silicate grains {e.g. Draine & 
Lee 1984) predicts two distinct grain temperatures. For particles with equal radii of 0.1 
/im, graphite and silicate grains are predicted to have temperatures of 18.8 K and 15.4 K, 
respectively. Taking into account the predicted abundance and absorption efficiency, the 
graphite grains produce about 60-75 % of the far-infrared emission. A two-composition 
model with the above mixture is a significantly worse fit to the far-infrared spectrum 
observed by FIRAS than a single component with free temperature. Further work is needed 
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to determine the appropriate temperature distribution of the grains. At this stage we can 
state that the FIRAS observations require that the range of temperatures that contribute 
to the far-infrared emission is very narrow, and does not allow distinct compositional 
components with comparable abundances. 



4.2. Very Cold Component (4-7 K) 

The emission that we have characterized as a modified blackbody is a new component 
of the interstellar spectrum. It would have been difficult to detect without the continuous 
spectral coverage of FIRAS, and it only becomes important at frequencies below 15 cm~^; 
see Figure ^. In the first paper presenting FIRAS observations of the Galaxy, Wright et al. 
found that the very cold emission was correlated with a galactic angular template based on 
the total intensity in the FIRAS spectra (over the cosmic background radiation). The very 
cold component was also found in the FIRAS data by Barnes (1994), who used a principal 
component analysis to separate the FIRAS data into 'eigenspectra' and 'eigenmaps'. We 
have now determined the angular distribution of the very cold component by modelling 
coadded spectra regions at high and low latitude. The salient results for constraining 
the nature of the very cold component are the following: (1) The temperature of the 
very cold component varies little with galactic longitude, and is not correlated with the 
temperature of the warm dust. (2) The very cold component is ubiquitous, and its optical 
depth is correlated (though not directly proportional to) the optical depth of the warm dust 
emission. We consider now some hypotheses for the nature of the very cold component. 



4-2.1. Dust shielded from the interstellar radiation field? 

A straightforward hypothesis for the origin of the very cold component is dust heated 
by an attenuated radiation field. For example, dust inside molecular clouds with extinction 
Ay > 1 is cooler because it is protected from ultraviolet photons. In order to explain the 
low temperature of the very cold component, the heating rate must be attenuated by a 
factor (T1/T2Y ~ 10'^. Calculations of the penetration of the interstellar radiation field 
into molecular clouds indicate that the heating rate is attenuated by a factor of 10^ at the 
center of a cloud with extinction Ay ~ 20 for silicate grains and Ay ~ 50 for graphite 
grains (Mathis et al. 1983, hereafter MMP). Similar results are obtained by somewhat more 
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detailed models that take into account the self-heating of the cloud by infrared emission 

from the outer layers (Bernard et al. 1992). The corresponding column density of gas is 
Nh ~ 4 — 10 X 10^^ cm~^ (Bohhn, Savage, & Drake 1978). By the standards of the diffuse 
interstellar medium, this is a large column density, so clouds containing shielded dust 
should be well known. 

We note that a much lower estimate for the column density needed to shield dust to 
temperatures corresponding to the very cold component was obtained by Barnes (1994), 
who concluded that the very cold component could be explained by dust shielded by the 
outer layers of molecular clouds with Ay = 3 mag. Barnes scaled the intensity of starlight 
in clouds by the exponential factor, e~'^'^^ , where tuv is the far-ultraviolet optical depth. 
This estimate neglects the ability of visible photons, which contain more than half of the 
energy density of the interstellar radiation field, to penetrate much deeper than UV photons 
(cf. Mathis et al. 1981). Further, scattering allows even ultraviolet radiation to penetrate 
much deeper than the exponential factor predicts; for example, starlight in the center of a 
cloud with total extinction Ajjy = 20 is predicted to be ~ 50 times brighter when scattering 
(with uj = g = 0.5) is included (Flannery, Roberge, & Rybicki 1980). Very large column 
densities are needed to shield dust from the interstellar radiation field. 

The very cold component is observed to be ubiquitous; in particular, it is present in the 
high-latitude spectra. In order for shielded dust to explain the very cold component, there 
must be clouds with high extinction distributed throughout the sky. This possibility can be 
excluded, and we can show that their filling factor must be small on angular scales down to 
5'. Because we can see external galaxies, it is clear that dark clouds cannot fill much of the 
sky on scales larger of a few degrees. The Shane & Wirtanen (1967) galaxy counts within 
1 deg^ bins are relatively smooth, and are anticorrelated with the 21-cm line brightness 
(Heiles 1976). The bin-to-bin dispersion in the galaxy counts (~ 25% at \b\ > 20°) is due 
primarily to the 'intrinsic mottling in the surface density of the distribution of galaxies' 
(Burstein & Heiles 1978). Thus on very large scales, there is no evidence for optically thick 
clouds that could shield dust from the interstellar radiation field. 

On intermediate scales, regions containing shielded dust should be observable as dark 
clouds against background starlight. The Palomar Observatory Sky Survey has been 
surveyed for all regions of noticeable extinction (Lynds 1962; Magnani, Blitz, & Mundy 
1985), and the regions of high-latitude extinction were found to be molecular clouds traced 
by CO rotational lines. However, the total extinctions of these clouds are generally small 
{{Av) — 1 for individual clouds; de Vries & Magnani 1986). Further, the covering fraction 
of the CO-emitting portion of high-latitude molecular clouds is small (~ 5 x 10~^; Magnani, 
Lada, & Blitz 1986; but see also Heithausen et al. 1993) so that their column density 
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averaged over large regions is only of order (Nh) ~ 10^^ cm~^. We may relate our derived 

optical depths at z/q = 30 cm~^ (333/im) to column densities using Tq = 1.3 x 10~^^A^j^ (cm~^) 
(based on 250/im observations; Hildebrand 1983). The observed optical depth of the cold 
component at |6| > 30° corresponds to a column density of (Nh) ~ 5 x 10^'^ cm~^. This is 
some two orders of magnitude greater than the amount of molecular gas at high latitudes. 

On angular scales down to the 5' resolution of the Infrared Astronomical Satellite 100//m 
sky survey, regions of shielded dust should be bright infrared emitters. Shielded regions 
must reemit all of the energy absorbed from the interstellar radiation field; for almost 
any solid material, this emission occurs at mid- to far-infrared wavelengths. The surface 
brightness of the optically thick regions, integrated over infrared wavelengths, must equal 
the surface brightness of the interstellar radiation field (cf. Keene 1981; MMP). The power 
in the very cold emission is observed to be negligible, so the energy would be emitted by 
the warm component. The range of temperatures observed for the warm component is not 
large, so the surface brightness is traced (within a factor of 2) by the lOO/im emission. The 
surface brightness of the ISRF is AnJ = 2.17 x 10~^ erg cm^^ s^^ (MMP), which corresponds 
to a specific surface brightness at lOO/im for T = 18 K and a = 2 of /,^(100//m) = 43 MJy 
sr^^. It is straightforward to see all regions with /,^(100/im) > 1 MJy sr^-*^ on the IRAS sky 
maps. The bright IRAS-detected clouds correspond to known dark clouds and high-latitude 
molecular clouds, with some exceptions (cf. Blitz, Bazell, & Desert 1990). As discussed 
above, the amount of material in high-latitude molecular clouds is far too low to explain 
the observed brightness of the very cold component of the FIRAS spectra. 

On small angular scales, the shielded regions would be visible as extinction on 
high-resolution optical images and would be point sources to IRAS. Optical images of 
high-latitude clouds reveal small-scale structure that does not always correlate with the 
infrared emission (Paley et al. 1991, Guhathakurta 1994). The fact that these regions 
are seen as bright on optical images suggests that the light is a combination of scattered 
starlight and luminescence. The lack of a 100/im counterpart precludes these regions being 
optically thick to the ISRF, so we agree with previous suggestions that the anomalous 
optical brightness is due either to an admixture of large grains or to radiative transfer 
[e.g. shadowing) effects. The existence of small, shielded regions can be ruled out on 
thermodynamical grounds. Such small, shielded regions would have substantial overpressure 
and would rapidly evaporate without a massive central object to hold them together: the 
pressure of shielded regions smaller than 5' within 100 pc is at least two orders of magnitude 
larger than the pressure in the diffuse ISM (cf. Reach, Heiles, & Koo 1993). 
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4-2.2. Very small grains? 



Very small grains (< 0.02//m) undergo temperature fluctuations when exposed to an 
ultraviolet radiation field, because the mean interval between successive ultraviolet photons 
is longer than the cooling time. Immediately after an ultraviolet photon is absorbed, very 
small grains cool by radiating at frequencies far higher than the larger grains. In the 
time interval between ultraviolet photon impacts, the temperature of the particle is very 
low, but at least the cosmic microwave background temperature (2.7 K). Nonequilibrium 
emission has been invoked to explain the mid-infrared emission of interstellar dust (Draine 
& Anderson 1985; Weiland et al. 1986; Desert et al. 1990). For example, a 0.001 //m silicate 
grain is predicted to reach maximum temperatures ~ 200 K, but to spend the bulk of the 
time at temperatures < 10 K (Draine & Anderson 1985). The emission from very small 
grains was modeled by Desert et al. (1990) from the near-infrared to 1 mm wavelength. In 
their model some submillimeter emission is predicted from the small grains, and it accounts 
for 12% of the emission at 800 fira. Based on our two-component fits to the spectra of 
the outer galactic plane (270° > / > 90°), the very cold component contains 1/3 of the 
total emission at 800 /im. If the models have underpredicted the 800 /xm emission of very 
small grains by a factor of ~ 3, for example because of our poor understanding of their 
low-temperature heat capacity, then very small grains could produce the cold component 
observed in the FIRAS spectra. However, stochastically heated particles producing the 
800 /im emission will produce excess emission at wavelengths shorter than 60 /im, and will 
produce mid- infrared colors inconsistent with IRAS and DIRBE observations. If there were 
a heating mechanism capable of maintaining very small grains at a minimum temperature 
~ 6 K, where they spend the bulk of their time, then they would be able to explain the 
very cold component; such a mechanism is not presently known. 



4.2.3. Grains with unusual optical properties? 



A population of cold dust grains, either well mixed with the warm grains or in 
widely-distributed clouds, could explain the very cold component. Balancing heating and 
cooling, the grain temperature scales as — T^^b. {Qvis/QFm)'JisRF, where Qvis is the 
absorption efficiency averaged over the interstellar radiation field, Qfir is the absorption 
efficiency averaged over the grain's emitted spectrum, Jisrf is the strength of the interstellar 
radiation field, and we have assumed a — 2. Particles with Qvis/QriR a factor of 10~^ lower 
than that of the grains responsible for the warm (T ~ 18 K) dust emission could explain 
the very cold component. 
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Practal grains have enhanced efficiency for submiUimeter emission, and therefore they 
can cool to very low temperatures (Wright 1993). A population of fractals may exist 
such that, when averaged over their size and shape distributions, the peak emission has 
characteristic temperatures corresponding to the observed cold component. Because of their 
high efficiency at emitting in the far-infrared, even relatively large fractal grains undergo 
temperature fluctuations (Bazell & Dwek 1990, Wright 1993). A population of fractal grains 
that spend the bulk of the time at very low temperature (not much above 2.7 K) and pulse 
to 5-10 K after absorbing a photon from the radiation field could produce the observed cold 
component of the FIRAS spectra. It is very important that the absorption cross-section per 
unit mass of fractal grains is substantially higher than for compact spherical grains. Thus 
the mass of fractal dust required to produce the cold component is not large. In order to 
explain the correlation of the warm and cold components, the fractal grains would have to 
be spatially correlated with the grains responsible for the warm component. 



4 .2. 4- Very large particles? 



Particles sufficiently large that they can cool via photons with wavelengths much 
smaller than their radii reach an equilibrium temperature in the interstellar radiation field 
such that T — {J/gsb + T^brY^'^ — 2-9 K, where J is the angular- averaged interstellar 
radiation field and ass is the Stefan-Boltzmann constant. Smaller particles emit with 
decreased efficiency in the far-infrared, and thus are warmer. Grains of astronomical silicate 
and graphite (Draine & Lee 1984) equilibrate at 6 K if they have radii of order 100 /im. 
Similar particles {e.g. 30;um amorphous C spheres) have been proposed before, in order 
to explain excess submiUimeter absorption in H II regions (Rowan- Robinson 1992). A 
problem with the hypothesis that large grains produce the very cold component is that the 
temperature of the cold component does not change much in the inner Galaxy, despite a 
substantially higher heating rate there. 

The mass of large grains required to explain the observed submiUimeter emission is 
not large. At a wavelength of 800 /^m, the absorption efficiency of 100 /im spheres is 0.2 
for graphite or 0.4 for silicate (Draine & Lee 1984). The mass of dust required to explain 
the very cold component is about 0.1% of the mass of interstellar gas. If the large grains 
contain G or (MgSiFe)04 groups, their abundance must be less than ~ 0.0073 relative to 
H by mass (based on the abundances in Anders & Grevesse 1989); this upper limit is a 
factor of 7 times larger than the inferred mass of large grains. Based on the strength of the 
10/xm bands seen in absorption toward bright infrared sources, nearly all of the cosmic Si 
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abundance is locked in grains smaller than l//m; similarly, infrared absorption bands and 
the 2200 A bump require ~ 60% of C to be locked in small grains (Tielens & AUamandola 
1987). This still leaves enough room for the large grains to be composed of either silicates 
or carbonaceous material. 

Size distributions of the form dn/da oc a~"*, up to a maximum size a^axi have been 
found to match the wavelength dependence of the extinction of starlight, where m = 3.5 
with flmax- = 0.25 yum (Mathis, Rumpl, & Nordsieck 1977; Draine & Anderson 1985). If this 
distribution were to extend from amax to 100 fim, there would be 50 times more particles 
than are needed to explain the very cold component of the far-infrared spectrum; therefore 
the slope of the size distribution from Qmax to 100 /xm would have to steepen to m ~ 4.2. 
We note that there are other limits to the number of very large particles. For example, 
inversion of the wavelength-dependence of polarization into a size distribution shows that 
the size distribution must steepen at sizes larger than ~ 0.7 iim (Kim & Martin 1994). Any 
model for the large particles would have to match the extinction and polarization curves, 
while remaining consistent with cosmic abundances. 



4-2.5. An emissivity enhancement? 



The cold component we observe may not be cold dust at all, but rather a spectral 
feature in the long-wavelength emissivity of the particles that produce the bulk of the 
far- infrared emission ('warm dust'). Such an emissivity feature could arise from resonances 
in the grains, impurities, or molecules in the mantles. The fact that the mixing ratio of the 
cold and warm components is nearly the same at every position motivates this hypothesis. 
In order to test this hypothesis, we divided the observed spectra by the warm component 
of the spectral models with z/^ emissivity at low frequency. This ratio, the emissivity 
enhancement, should reveal the shape of the spectral feature. The observed spectrum 
deviates significantly from the warm component at frequencies below 25 cm~^, and the 
excess at 12 cm~^ is about 30%. The enhancement increases toward lower frequencies as 
at least to 10 cm~^. 

If the excess submillimeter emissivity were due to low-frequency resonances in the 
grains, then it would be expected to eventually disappear at frequencies below the 
lowest resonant mode. As a check, we used the DMR observations to check that the 
two-component modified blackbody spectrum does not exceed the observed sky brightness 
at frequencies below those to which FIRAS was sensitive. After subtraction of the dipole 
component of the cosmic microwave background, the correlation slopes of surface brightness 
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with csc|6| (for |6| > 15°) were determined in the three DMR wavebands; the results are 
shown in Table ^ The FIRAS data in each frequency bin treated similarly, yielding the 
spectrum for a cosecant Galaxy model. This FIRAS spectrum was fitted by two modified 
blackbodies of the form in Equation ^, with emissivity index a = 2. The warm component 
has a temperature Ti = 17.72 ± 0.14 K and optical depth n = (1.74 ± 0.06) x 10-^ 
and the very cold component has a temperature T2 = 6.75 ± 0.23 K and optical depth 
T2 = (1.23 ± 0.13) X 10~^. The extrapolation of a one-component fit and a two-component 
FIRAS fit down to the DMR frequencies is listed in Table ^. Neither dust model exceeds 
the Galactic emission observed by DMR. The DMR observations cannot strongly constrain 
the low-frequency dust spectrum, because synchrotron emission, free-free emission, and 
cosmic microwave background anisotropies all contribute to the signal. But the fact that 
the two-component model does not over-predict the low-frequency emission means that it 
is possible that the excess dust emissivity extends well into the microwave band. More 
sensitive microwave observations of the diffuse galactic emission, and a more detailed study 
combining the radio, microwave, far-infrared, and optical (cf. Reynolds 1992) constraints 
on Galactic and cosmic sky brightnesses, are needed to determine the diffuse microwave 
emission of interstellar dust. 



4.3. Evidence for other components 

The spectra of the inner Galaxy are poorly fit even by two modified blackbodies. There 
is a residual emission that peaks around 45 cm~^ (220 /im), suggesting a temperature 
around 14 K. This emission is present in spectra within 40° of the galactic center, so it 
emanates predominantly from within 5.5 kpc of the galactic center. Comparing the variation 
of the brightness of this component to that of the CO(l-O) line from the Goddard-Columbia 
survey (Dame et al. 1987) reveals a rough correspondence, suggesting that this residual 
emission may be due to dust associated with molecular clouds. The temperature of this 
component is similar to that obtained from far-infrared spectra of two Bok globules with no 
local heating sources (Keene 1981). Thus dust shielded from the interstellar radiation field 
is present in the FIRAS spectra after all, at a temperature intermediate between that of the 
'warm' dust and that of the 'cold' dust. The heating rate of the intermediate-temperature 
dust is decreased from that of the warm dust by a factor ~ 5, suggesting extinctions Ay ~ 2 
between the intermediate-temperature dust and the interstellar radiation field. By fitting 
a three-component spectral model to the data, we found the typical optical depth of the 
intermediate-temperature component is ~ 1 x 10^^ at 30 cm^^, corresponding to a column 
density ~ 1 x 10^^ cm~^. The column density (inferred from the optical depth) and the 
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extinction (inferred from the temperature) are in good agreement with observations of 
A^(H) + 2A^(H2) and Ay for interstellar gas and dust (Bohlin, Savage, k Drake 1978). This 
component is very weak, contributing only ~ 2% of the emission at 45 cm~^ and ~ 60 times 
less power thatn the warm component. 



4.4. Compcirison to other galetxies 

Submillimeter and far-infrared observations of other galaxies using ground-based 
and airborne telescopes are challenging with current technology, and continuous spectral 
coverage of the dust emission comparable to that obtained by FIRAS for our Galaxy has 
never been achieved. But a comparison to observations of external galaxies is important 
for two reasons. First, the entire far-infrared emission of external galaxies can be measured 
unambiguously — without confusion from other emission sources such as zodiacal light and 
background radiation. And second, observations of other galaxies allow us to determine 
whether the far-infrared spectral shape is a generic property of cosmic dust or is peculiar 
to individual galaxies. Most galaxies observed in the submillimeter range are significantly 
warmer than our Galaxy. Eales, Wynn- Williams, & Duncan (1989) observed a sample of 10 
galaxies (with large 100/xm fluxes) at 350/xm, with some observations at 450, 800, & 1100 
/im as well. Their beamsize corresponded to 3 kpc for the nearest and 30 kpc for the most 
distant galaxies observed. The FIRAS spectra of the galactic center and inner Galaxy are 
substantially cooler than the average spectrum of Eales et al. : the submillimeter fluxes 
are relatively higher and the 60/um flux is lower. In another survey, Clements, Andreani, & 
Chase (1993) reported observations of 5 galaxies at 450, 800, & 1100 /xm, with beamsizes 
ranging from 1.4-12 kpc. For the 5 galaxies, the temperature of a single-component flt 
ranged from 28-35 K, compared to the range 16-23 K for spectra of our Galaxy. Recently, 
Franceschini & Andreani (1995) reported on a 1.25 mm survey of a complete sample of 
southern galaxies, with a selection criterion based on the IRAS 60//m flux. Their average 
spectrum is also 'warmer' than that of the Milky Way, and they found no evidence for 
substantial emission from cold dust. Evidently the galaxies selected in these surveys are 
warmer than our Galaxy because they were selected based on their brightness at 60-100 yum: 
Using the IRAS surface brightness to select galaxies apparently selects the warmest galaxies. 
Based on the observed temperatures, it is likely that the average radiation field within these 
galaxies is nearly an order of magnitude stronger than that in our own Galaxy. 

One other study has shown that the far-infrared spectrum of the Milky Way is not 
atypical. Stark et al. (1989) observed three spirals in the Virgo cluster at 160/xm and 
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360/im with a resolution corresponding to 6.5 kpc. When we combine the FIRAS spectra 
within 30° of the galactic center, a region that should correspond to what Stark et al. would 
have observed if our Galaxy were at the distance of Virgo, we find that the spectrum of our 
Galaxy is entirely consistent with that of the Virgo spirals. 



4.5. Implications for Cold Gas as Dark Matter 

The lack of substantial power emitted by dust with temperatures in the 4-15 K range 
places strong constraints on the amount of very cold material in the interstellar medium. 
The discovery of CO absorption lines from clouds in the outer galaxy led to the suggestion 
that very cold molecular clouds, which would be very difficult to detect by CO emission, 
may be an important, massive component of the Galaxy (Lequeirx, Allen, & Guilloteau 
1993). The power of all starlight incident on such cold clouds would be reemitted within the 
FIRAS bandpass. In fact, there is little power emitted besides that contained in the warm 
component; the fraction of the total power contained in the cold component is shown in 
Figure ||. Because the cold component is detected in every direction, it cannot be attributed 
to outer galaxy molecular clouds; therefore the power shown in Figure |^ is an upper limit. 
Further, the cold component contains less than 0.05% of the power in the local interstellar 
radiation field, so that it cannot be due to optically thick clouds if it is produced locally. 

Any model for the Galactic interstellar medium including a substantial increase in the 
abundance of dark clouds must comply with the constraint that very little power is emitted 
by cold dust. This constraint has been neglected in some recent work {e.g. Lequeux et al. 
1993; Pfenniger, Combes, & Martinet 1994a). One possibility is that such clouds have a 
very low dust abundance. This would be consistent with the low gas temperature, because 
the photoelectric heating by dust grains would be absent. However, there is evidence for 
star formation in the far outer galaxy (De Geus et al. 1993), so there should be dust 
formation in stellar outfiows and supernovae. Another possibility is that the radiation field 
is a factor of ~ 10^ lower in the region of the cold clouds. This too is contradicted by the 
presence of star formation. Furthermore, power absorbed from the extragalactic radiation 
field would still be a factor of ~ 10 brighter than our upper limit. 



5. Conclusions 
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The continuous spectrum of galactic dust observed by the COBE/FIKAS can be 
characterized as a sum of three emission components: warm dust heated by the interstellar 
radiation field, intermediate-temperature dust heated by the attenuated radiation field in 
molecular clouds, and a very cold component. The very cold component is correlated with 
the warm dust, emits very little power, and its temperature does not increase in the inner 
Galaxy. The very cold component is unlikely to be due to dust that is cold because it is 
shielded from the interstellar radiation field. The very cold component may be a spectral 
feature in the warm dust constituent material, or may be due to fractal or needle-like grains 
or particles substantially larger than those responsible for the warm component. The lack 
of substantial power from cold dust limits the amount of cold gas in the outer galaxy and 
makes it unlikely that such gas could contribute an unseen, massive component of the 
Galaxy. 

The existence of enhanced submillimeter emission from interstellar dust, regardless of 
its origin, implies an increased contribution to the microwave sky brightness by the Galaxy 
than previously expected from extrapolations of 100 fim observations {e.g. by IRAS). The 
interstellar dust emission can be characterized by a two-temperature approximation with 
temperatures and optical depths as shown in Figure ^. The implications for searches for 
fluctuations in the microwave background radiation are not yet clear. If the cold component 
is due to enhanced emissivity of the warm dust constituent material, then its spatial 
fluctuations will be comparable to those of an IRAS lOO/xm map. However, if the cold 
component is due to a separate population of dust grains in clouds with a substantially 
different morphology than previously considered, the fluctuations in the microwave sky 
brightness due to interstellar dust could be notably enhanced. 

WTR thanks M. G. Hauser and S. H. Moseley for insightful discussions. 
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Table 1: Optical depths and temperatures in the galactic plane 
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Table 2: Optical depths and temperatures at high latitude 
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20.7 ±0.5 


1.36 ±0.06 


2.37 


0.95 


135 </ < 180, 10 < 6 < 30 


0.60 ±0.02 


19.3 ±0.2 


1.51 ±0.03 


1.88 


> .99 


135 < Z < 180, -10 > 6 > -30 


0.88 ±0.01 


19.1 ±0.1 


1.60 ±0.02 


2.85 


> .99 


180 < Z < 225, 10 < 6 < 30 


0.20 ±0.01 


22.0 ±0.6 


1.12 ±0.06 


1.78 


> .99 


180 < Z < 225, -10 > 6 > -30 


0.79 ±0.01 


20.8 ±0.1 


1.58 ±0.02 


2.68 


> .99 
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Table |[ (continued) 



(c) Two- component fit (a = 2) 



region 


lOVi 


Ti (K) 


10^2 


T2 (K) 


xVdof 


P{> Fac) 


P{> Ffec) 


6 > 60 


0.13 ±0.01 


17.5 ±0.3 


1.2 ±0.2 


7.1 ±0.3 


4.50 


> .99 


0.50 


< / < 90,30 < 6 < 60 


0.19 ±0.01 


17.8 ±0.2 


1.2 ±0.1 


7.3 ±0.3 


3.02 


> .99 


0.48 


< / < 90, -30 > 6 > -60 


0.25 ±0.02 


17.5 ±0.2 


1.4 ±0.4 


6.7 ±0.6 


2.31 


0.96 


0.56 


135 </ < 180,10 < 6 < 30 


0.83 ±0.02 


16.8 ±0.1 


4.9 ± 1.0 


5.7 ±0.3 


1.66 


> .99 


0.76 


135 </ < 180,-10 > 6 > -30 


1.16 ±0.01 


17.0 ±0.1 


9.5 ± 1.3 


5.1 ±0.2 


1.47 


> .99 


> .99 


180 </ < 225, 10 < 6 < 30 


0.34 ±0.02 


17.0 ±0.2 


3.3 ±0.8 


6.1 ±0.4 


1.72 


> .99 


0.58 


180 </ < 225,-10 > 6 > -30 


1.04 ±0.01 


18.3 ±0.1 


7.9 ± 1.0 


5.5 ±0.2 


1.70 


> .99 


> .99 
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Table 3: Comparison of dust models to DMR observations 

FIRAS dust models 

frequency DMR csc|6| slope one-component two-component 

(GHz) (cm-i) (kJy sr"!) (kJy sr'^) (kJy sr'^) 

31 To 1.5 ±0.2 Om 004 

53 1.8 1.3 ±0.2 0.1 0.3 

90 3.0 2.5 ±0.8 0.9 2.2 
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Fig. 1.— The FIRAS spectrum of the region (85° < I < 95°, 25° < 6 < 35°) is shown together 
with its decomposition into cosmic background radiation (CMBR), Galactic emission, and 
zodiacal light. 

Fig. 2. — (a) Spectrum of interstellar emission in the galactic plane toward longitude / = 45°. 
Other than the bright spectral lines arc due to and N"^, the spectrum is dominated by 
emission from warm dust, which peaks around 65 cm~^ (150/im). (b) Residuals after a 
single-component modified blackbody with a — 2 was subtracted. The residual intensity 
at each frequency was divided by the uncertainty at that frequency. The excess emission 
at 7-20 cm~^ (1400-500/im) is evident, (c) Residuals after a single component with the 
best-fitting emissivity index was subtracted. The fit has improved with the extra degree 
of freedom, but the true shape of the spectrum at low frequencies has not been accurately 
fitted, (d) Residuals after a two-component model was subtracted. The very cold emission 
is now adquately fit. The spectral lines of CO(4-3) and CI can be clearly seen. There is 
evidence for excess emission in the 40-55 cm~^ (250-180/im) for the inner galaxy spectra. 

Fig. 3. — (a) Spectrum of interstellar emission in the region with —30 > b > —60° and 
180° > I > 90°. Only a weak C+ spectral fine and the warm dust continuum are evident. 
(b) Residuals after a single-component modified blackbody with a = 2 was subtracted. 
The excess emission at 7-20 cm^^ (1400-500;um) is evident, (c) Residuals after a single 
component with the best-fitting emissivity index {a = 1.4) was subtracted, (d) Residuals 
after a two-component model was subtracted. It is not possible to clearly distinguish whether 
panel (c) or (d) are better fits to the data. 

Fig. 4. — Longitude profiles of the optical depth and temperature of the two-component, 
four-parameter fits to the spectra of the galactic plane. 

Fig. 5. — Maps showing the variation of the temperature and optical depth of the 'warm' 
and 'cold' components. The maps are in galactic coordinates, with the galactic center in the 
middle. The grid lines delineate the regions used for spectral fitting, and the values listed 
within each region are the fitted parameters and their statistical uncertainties. Across the 
galactic plane, the longitude variation from the previous Figure is repeated for comparison. 
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Fig. 6. — Optical depth of the cold component versus optical depth of the warm component 
for all spectra. The two components are well correlated, which strongly suggests that the 
cold component is galactic emission. Because the galactic column density varies strongly 
with latitude, the points in the lower- left-hand corner are all at high galactic latitudes. The 
cluster of points in the upper-right-hand corner are the galactic plane spectra; the outer 
galactic spectra are at the lower left of this cluster and the galactic center is at the top right. 
The line representing T2/T1 — 7 is shown for comparison. 



Fig. 7. — The FIRAS galactic spectrum of the region (/ = 45°, \b\ < 3°) is shown together 
with its decomposition into warm and very cold components. 



Fig. 8. — Fraction of the total power contained in the 'cold' component, as a function of 
galactic longitude. 
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